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Abstract. The decay of an excited atom in the presence of a medium that 
both scatters and absorbs radiation is studied with the help of a quantum- 
electrodynamical model. The medium is represented by a half space filled with 
a randomly distributed set of non-overlapping spheres, which consist of a linear 
absorptive dielectric material. The absorption effects are described by means of 
a quantized damped-polariton theory. It is found that the effective susceptibility 
of the bulk does not fully account for the medium-induced change in the atomic 
decay rate. In fact, surface effects contribute to the modification of the decay 
properties as well. The interplay of scattering and absorption in the total decay 
rate is discussed. 
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1. Introduction 

The spontaneous-emission rate of an excited atom can be altered by the atomic 
environment, as has been pointed out long ago [1]. For an atom embedded in a uniform 
linear non-absorptive dielectric of infinite extent the change in the emission rate has 
been obtained from quantum electrodynamics [H |3] . For dense media local-field effects 
have to be considered as well [H [5]. If the medium can absorb the emitted photons, 
the analysis gets more complicated, since the loss mechanism has to be treated in a 
quantum-mechanical context [B] - |14) . 

Local- field effects do not play a role if a different geometry is considered, with the 
atom situated outside a medium of finite (or semi-infinite) extent. A well-known case 
is that of an atom in front of a medium that fills a half-space. For this configuration 
the decay rate depends on the distance between the atom and the medium [15] - [22] . 

In all treatments mentioned so far the medium is structureless on the scale of the 
wavelength corresponding to the atomic transition. As a consequence, the medium 
properties are fully described by a susceptibility, which does not vary appreciably 
on the scale of the wavelength. The picture changes if the structure of the medium 
cannot be neglected, since scattering may occur then as well. The total extinction in 
such a medium is driven by both absorption and scattering. In practice, extinction 
by scattering in material media is quite common. Impurities and defects both lead to 
scattering effects that are difficult to avoid. 

To study the interplay between the two types of extinction that may modify 
atomic decay processes in the presence of a material medium, it is useful to analyse a 
model in which both of these features occur simultaneously. The model that we shall 
adopt in the following is that of a medium consisting of non-overlapping spheres that 
are made of an absorptive material. The spheres, which may move freely through the 
system, are distributed randomly with a uniform average density. In a recent paper [2^ 
a similar model with a collection of spherical scatterers consisting of non-absorptive 
material has been studied. 

In order to describe the absorptive dielectric material of the spheres we shall use 
the quantum-mechanical damped-polariton model. The central quantity in this model 
is a space-dependent polarization density, which is coupled to the electromagnetic field 
and to a bath of harmonic oscillators accounting for absorption. The Hamiltonian of 
the damped-polariton model can be diagonalized exactly, as has been shown both for 
the case of a uniform dielectric [24 and for a dielectric with arbitrary inhomogeneities 
[25]. 

To arrive at analytical results for the decay rate in the presence of a medium with 
extinction due to both absorption and scattering we will adopt several approximations. 
The density of the spherical scatterers will supposed to be low, so that the medium 
is dilute. Furthermore, the size of the spheres will be taken to be small as compared 
to the atomic wavelength. Finally, the distance from the excited atom to the medium 
will be chosen to be large compared to the wavelength. 

The paper is organized as follows. In section 2 the properties of the model and 
its diagonalization will be summarized. In section 3 the decay rate of an excited 
atom in the presence of an arbitrarily inhomogeneous damped-polariton dielectric 
will be derived from the basic Hamiltonian. The decay rate is determined by the 
electromagnetic Green function, which enters the description via a specific coefficient 
in the diagonalization matrix. Since the medium consists of randomly distributed 
spheres the Green function has to be averaged over all configurations in order to 
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obtain the physical decay rate. This averaging procedure wih be discussed in section 
4 and 5. As it turns out, the averaging process for a medium that fills a finite region 
of space (or a half-space) should be carried out carefully, since the boundaries give rise 
to specific surface contributions. Once these surface effects have been evaluated for 
the specific case of a medium filling a half-space, we can obtain the average decay rate 
of an excited atom in the presence of such a medium. The change in the atomic decay 
rate as a function of the distance between atom and medium will be determined in 
section 6, and the interplay of absorption and scattering processes will become clear. 
Some of the technical details of the derivation are given in two appendices. 



2. Field quantization in the presence of an inhomogeneous absorbing 
dielectric medium 

In the damped-polariton model the dielectric medium is described by a polarization 
density, which interacts with the electromagnetic field according to the standard 
minimal-coupling scheme. To account for absorption effects a bath of harmonic 
oscillators with a continuous range of eigenfrequencies is coupled to the polarization 
density in a bilinear way. The Hamiltonian of the model is [24j [25] 

f r 1 1 1 1 f°° 

H,^J,r [-tf + -(V X A)^ + -P^ + lp.iX^ + - 1 a.Ql 
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(2.1) 



The transverse part of the electromagnetic field is determined by the vector 
potential A(r), for which the Coulomb gauge is adopted. Its conjugate canonical 
momentum is n(r). The linear dielectric, with a space-dependent density p(r), is 
described by the harmonic displacement variable X(r) and its canonical momentum 
P(r), with the associated eigenfrequency a;o(r). The electromagnetic field is coupled 
to the dielectric variables in the usual way. In terms of the polarization density 
— Q!(r)X(r), with a space-dependent coupling parameter a{r) > 0, the minimal- 
coupling scheme leads to an electrostatic contribution involving [Q!(r)X(r)]i and to 
a bilinear interaction term with A(r) • P(r). The longitudinal part of a vector (or 
a tensor) is obtained by a convolution with the longitudinal delta function ^^(r) = 
— W(47rr)^^. Finally, damping is introduced in the system by a continuum bath of 
harmonic oscillators with canonical variables (r) , Q,^ (r) and with eigenfrequencies 
uj. These bath oscillators are coupled to X(r) with a strength v^{r) > 0. 

The canonical variables obey the standard commutation relations: 

[n(r),A(r')] = -indrir-r') [P(r),X(r')] = -ih\6{r-r') 

[Q^{r),Y^,{r')] ^ -ihS{LJ - oj')\d{r - r') (2.2) 

while all other commutators of the canonical variables vanish. Here I is the three- 
dimensional unit tensor, while ^T(r) = I S{r) — ^^(r) is the transverse delta function. 
The electric field operator E(r) is the sum of a transverse part depending on 11 (r) 
and a longitudinal part that is proportional to the polarization density: 

E(r) = -ln(r) + lKr)X(r)U. (2.3) 

£0 EO 
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The Hamiltonian is quadratic in the canonical variables, and can be diagonalized 
explicitly pS] : 



Hd = JdrJ^ dw/ia;Ct(r,a;) • C 



(r,a.) (2.4) 



where we omit a zero-point-energy term. The operators C(r,a;) are annihilation 
operators, which (together with the associated creation operators) satisfy the 
commutation relations: 

[C(r,a;),C^(r',c.')] ='5(c^-c.')l<5(r-r') [C(r, ..), C(r', c')] = . (2.5) 

Each canonical operator can be written as a linear combination of the annihilation 
and creation operators. For instance, one has 

A(r)=y"dr'^ dujfA{r,r',uj)-C{r',uj)+h.c. (2.6) 

E(r)=y"dr'^ dcjfi;(r,r',w)-C(r',Lj)-Kh.c. (2.7) 

with the tensorial coefficients f^i and f^;. Similar expressions can be written for the 
other canonical variables. Since the vector potential is transverse, the coefficient (a is 
transverse in r. 

In order to derive explicit expressions for the coefficients one may use a method 
due to Fano [IHI- It amounts to evaluating the commmutator of C(r, w) with the 
Hamiltonian Hd in two different ways. On the one hand this commutator follows 
from (|2.4p and (|2.5p as [C{r,uj), Hj] —hLuC{r,uj), and on the other hand it may be 
evaluated by first writing C(r,w) as a linear combination of the canonical variables, 
subsequently inserting (|2.1|) and finally employing p.2p . Upon solving the linear 
equations that follow by comparing the results of these two approaches, one arrives at 
explicit expressions for the coefficients in terms of the tensorial Green function G of 
the system [25]. The latter is defined as the solution of the standard equation 

^2 

^Vx [V xG{r,r\uj + iO)] + — [1 + x(r, w + iO)] G(r, r', + iO) = I 5(r - r') (2.8) 

with a; -|- iO in the upper half of the complex plane and infinitesimally close to the 
positive real axis. In the course of the diagonalization process the frequency- and 
position-dependent susceptibility x is found as 

-1 



a2 



x(r, i:j-|-iO) 

£oP 



2 2 



do; ■ 



.(2.9) 



" P^Jo (c.-t-iO)2-C.'2 

The tensorial Green function satisfies the reciprocity relation 

G(r,r',a; + iO) = G(r',r,a;-|-iO) (2.10) 

where the tilde denotes the tensor transpose. In terms of the above tensorial Green 
function, coefficients (j2.6p and (|2.7|) are given as 

^ / /X /Sim Y(r',a; + «0)\ ^ , , 

ffi r,r',cj^ -^— '-] G r,r',w + iO 2.11 

U[v,v',Lo) = --[iE[v,v',u)]T. (2.12) 

LO 

These expressions follow from the results presented in !25; . When they are substituted 
in (12. 6p and (|2.7p . they lead to expressions for the vector potential and the electric 
field that agree with those postulated in a phenomenological quantization scheme |11) , 

IIS, in. 
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3. Decay of an excited atom in the presence of an inhomogeneous 
absorbing dielectric 

When a neutral atom at a fixed position is present as well, the total Hamiltonian H 
of the system is given by the sum H^i + Ha + Hi of the damped-polariton Hamiltonian 
(|2.ip . the atomic Hamiltonian 



" Z^2m^^87reo|r.-r,| 2^ ineoh - ra\ ^ ' 

(with Z the atomic number, r^ the fixed position of the nucleus and r^ , the positions 
and momenta of the electrons) and the interaction Hamiltonian Hi , which follows from 
the usual minimal-coupling scheme as 

H,^ j Av[pa[v) ifiir) - Ja(r) • A(r)] . (3.2) 

Here pa{r) = eJ^iiH^ - ^a) - S{r - r^)] and Ja(r) = -^e X^JPi/™, '5(r - ^i)} 
are the local atomic charge and current densities, with curly brackets denoting the 
anticommutator. Furthermore, <p(r) is the scalar potential due to the polarization of 
the dielectric. Its gradient is given by 

V^(r) = -- [a{r) X(r)]L = -E(r) - - n(r) = -E(r) - A(r) (3.3) 

with the time derivative given as A = {i/h)[H, A]. In writing p.2p we have assumed 
that local-field effects are negligible. 

We assume that at the initial time t — 0, the atom is prepared in an excited 
state |e), while the dielectric medium (including the bath) and the electromagnetic 
field are in the ground state |0) of the Hamiltonian (|2.1[) or (12. 4p . i.e. in the state that 
is annihilated by all operators C(r,a;). The atom will decay to its ground state \g) 
with a time-dependent decay rate T{t). This rate follows from perturbation theory in 
leading order as 

= Tdt'eC/'^'^^/"^')*' mnrn' + cc. (3.4) 

n / "'0 

with \i) and |/) the initial and final states of the total system, with energies Ei and 
Ef, respectively. 

Upon taking the matrix element of p.2p . using charge conservation in the form 
V- jQ(r) — — (i//i)[iJa, yOa(r)], carrying out a partial integration and substituting (|3.3p 
with (I2.6l) - (|2.7p . we may rewrite the time-dependent decay rate as 

T{t) = ^5— / dt' / dr / dr' / dr" / dcj e'^'^"'^'')*' 
H ujI Jo J J J Jo 

(e|J,(r')|.g) • [(c., - c.)f^(r',r,L.) - if£(r',r,c.)] • 

(w,-u;)f^(r",r,t^)-|-if;(r",r,t^)] • (g|J„(r")|e) + c.c. (3.5) 

where Huja is the difference between the energies of the excited and the ground state 
of the atom. 

For large values of t the decay rate becomes independent of time. Carrying out 
the integrals over t' and w one finds that the coefficient drops out. As a result we 
get 

r = ^|dr| dr'l dr" (e| J.(r')|g) • fi=;(r', r, w.) • f;(r", r, c.,) • (5I J4r")|e) .(3.6) 
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Inserting (|2.1ip we obtain 

r=^y' dr'l dr"(e|J,(r')|5>-F(r',r",..„).(g|JJr")|e) (3.7) 
with the abbreviation 

F(r',r",tj) = y"drG(r',r,w + iO) •G*(r",r,a; + iO)Imx(r,a; + iO). (3.8) 

With the help of the differential equation (|2.8|) and the reciprocity relation (|2.10|) one 
may rewrite F(r,r',a;) as — (c/w)^ Im G(r, r', w + iO), so that the decay rate gets the 
final form 

^J-j I dr'l dr"(e|J„(r')|5).ImG(r',r",c., + iO)-(5|Ja(r")|e). (3.9) 

In the electric-dipole approximation the matrix element (e|JQ(r)|(7) is replaced by its 
localized form (5(r-ra) / dr(e|Ja(r)|g) = i (5(r-ra) WQ(e|/x|.g), with /x = -e X]i(ri-ra) 
the electric dipole moment. In that approximation the decay rate reads 

2(^2 

r = {e\ti\g) ■ ImG(r„ r„, + iO) • (g|/x|e) . (3.10) 

This expression for the decay rate, which is valid for an excited atom in the presence 
of an absorptive dielectric with arbitrary inhomogeneities, can be obtained as well by 
invoking the fluctuation-dissipation theorem [9l [13] . The above derivation shows how 
it follows from the explicit diagonalization of the inhomogeneous damped-polariton 
model in a rigorous way. 



r = -- 



4. Medium of absorbing spherical scatterers 

Let us consider a medium of non-overlapping spheres of absorptive material. It 
is an inhomogeneous dielectric that may be described by Hamiltonian (j2.1[) . The 
susceptibility (j2.9|) has a constant value within the spheres and vanishes outside, so 
that it may be written as x(r, uj + iO) — + iO) / (r)- If the radius of the spheres is 
a and the centre of the sphere i is located at r^ (with jr^ — rj| > 2a for i ^ j so as to 
avoid overlap), the function /(r) equals 9[a—\r — Yi\), with the step function 9[x) 
equal to 1 for a; > and elsewhere. 

We are interested in the decay of an excited atom in the presence of such a medium 
of absorptive spheres. Since the spheres may move the effective decay rate follows from 
(|3.9P or ()3.10|) by averaging over the positions of the centers of the spheres. Hence, 
we have to find an expression for the (imaginary part of the) average Green function. 

The differential equation (|2.8p for the Green function G(r, r', w -I- iO) is equivalent 
to an integral equation that relates G to the vacuum Green function Gq: 

G(r, r', z) - Go(r, r', z) - ^ x{z) j dr" /(r") Go(r, r", z) • G(r", r', z) (4.1) 

with the frequency variable z = + iO. The vacuum Green function is the solution 
of (|2.8I) with x(r,w -t- iO) = 0. Iterating (14. ip we get a series of terms, which up to 
second order in the susceptibility reads 

G(r, r', z) = Go(r, r', z) - ^ x{^) j dr" /(r") Go(r, r", z) • Go(r", r', z) 

+ ^ [X{z)? j dr" j dr'" /(r") /(r'") Go(r, r", z) • Go(r", r'", z) • Go(r'", r', z) + . . . 

(4.2) 
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When averaging both sides of this equation over the positions of the centers of 
the spheres, we need expressions for the averages (/(r)) and (/(r) /(r')). When the 
centers of the spheres are uniformly distributed with the density n, these averages 
have the form 



(/(r)) =nvo 
{f{v)f{v'))=n 

dr" ( dr' 



dr" 6l(a - |r-r"|)6i(a- |r'-r"|) 

'){a ~\v~ r"|) e{a - |r' - r"'|) .g(r", r"') 



(4.3) 



(4.4) 



with = Atto^ /i being the volume of the spheres and ^(r, r') being the pair correlation 
function. If the spheres are dilutely distributed, the correlations may be neglected, so 
that g can be replaced by unity. Upon carrying out the geometrical integrals expression 
(|4.4p becomes 

(/(r) /(r')) ^n[v^~ i^c? |r - r'| + ^ ^ |r - rf ] 0(2 a - |r - r'|) + vl . (4.5) 

Inserting the above averages in the iterated integral equation (|4.2p . we get up to 
second order in the susceptibility: 

Go(r,r',z)-^ni;ox(^) ^ dv" ^ dr'" Go(r,r",z)- 



(G(r,r',z)) 



I S{r" - r" 

^2 



■x{z)nvoGo{r",r" 



x(^)c(|r"-r"'|)Go(r",r"',z) 



Go(r" 



with the abbreviation 



e{2a 



4a ■ 16a3 

The right-hand side is the iterated solution of the integral equation 



(4.6) 
(4.7) 



(G(r,r',z)) = Go(r,r', 



- ^ I dr" I dr'" Go(r, r", z) ■ Xe{r" , r'", z) ■ (G(r"', r', z)) 

(4.8) 

again up to second order in x{z)- The effective susceptibility tensor is given as 



Xe(r,r',z) ^nvoxiz) 



l<5(r-r')-^x(z)c(|r-r'|)Go(r,r',z) 



(4.9) 



It should be noted that the term proportional to [x{z)]'^ in (|4.6|) results upon 
iterating (j4.8p up to second order in Xe- 

The effective susceptibility (j4.9|) is non-local with a range equal to 2a. The Green 
functions in the integrand of the last term of (|4.8p do not change appreciably over that 
range, when a is small compared to c/uj (for z = oj + iO), and to |r — r"| and |r' — r"'\. 
The first of these conditions can easily be met for spheres that are sufficiently small. 
In fact, we shall use (|4.8I) for oj equal to the atomic transition frequency, so that c/cj 
equals the transition wavelength. In contrast, the last two conditions are fulfilled only 
when r" and r"' in (14. 8p are sufficiently far from the fixed positions r and r'. Since 
the integrations in ()4.8p extend over all parts of space accessible to the spheres, it is 
not obvious that the two conditions can be fulfilled. We will postpone a discussion of 
this point to the end of the section. 
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When all three conditions mentioned above are satisfied, one may replace Xe by 
its localized version. Quite generally the localized version of a function F{r) that is 
of short range and centred around the origin can be written as a series expansion of 
which the first few terms are 

F(r) = ^ (r) J dr' F(r') - [V (5(r)] • J dr' r' i^(r') + | [ VV S{r)] : J dr' r'r' F(r') + . . . 

(4.10) 

We can evaluate the first few moments of F(r) = c(|r|) Go(r,0,w + iO) for small values 
oi q = oja/c by employing the expression [371 HH] 

Go(r, 0, z) = -— (\-—) e'^''/'' + V— ( -i — + f I - 3—) e'^''/'' 

+^^(r)l (4.11) 

for the vacuum Green function (with z in the upper part of the complex plane). The 
principal- value sign denotes the exclusion of an infinitely small spherical volume in 
subsequent integrations. As a result we obtain the following for the localized form of 

Xe- 



^ ^ I 1 4 2 2i 3 

1 — xiz) a 

^ * 3 15 ^ 9 ^ 



I 6{r - 



-nvo [x{z)? ^ (I A - 3VV) 6{r - r') . (4.12) 
75 

Here we have used the integral identities / dQ.rirj/r^ — (47r/3)(5ij and 
/ dVLrirjrj,ri/r^ = {^tt /lb) [6ij 6m + ^ifc + SuSjk), with being the cartesian 
components of the position vector r and dfl being an element of solid angle in the 
direction of r. 

When the localized form (I4.12p is inserted in (14.81) and a partial integration with 
respect to r" is carried out, we can employ the identity 

(I A - 3 VV) • Go(r, r', z) = Go(r, r', z) (4.13) 

for r 7^ r', as follows from the differential equation (|2.8p (with x = 0) for the vacuum 
Green function. It should be noted that the second operator between the brackets in 
(|4.13p does not contribute for r 7^ r', as is obvious from the form of the differential 
equation. There is no need to discuss the form that (|4.13p may take for r = r', since 
the localized form of the effective susceptibility can be used only when the arguments 
of the Green function are sufficiently far apart, as we have seen above. 

When the identity (|4.13p is taken into account, the localized form (I4.12p of the 
effective susceptibility may be rewritten as 



nvQxiz) 



l<5(r-r') =Xe(^)l'5(r-r'). 

(4.14) 

This localized form of the susceptibility has to be inserted in the integral equation 
(j4.8|) . Its solution up to second order in xi^) and up to first order in n reads 

(G(r, r', z)) = Go(r, r', z) - J Xe{z) J dv" Go(r, r", z) ■ Go(r", r', z). (4.15) 
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In principle, this expression for the average Green function, with r = r' = 
and z — (jja + iO, could be used to evaluate the right-hand side of (|3.10p . When the 
medium of scattering and absorbing spheres is infinitely large, the atom is necessarily 
embedded in the medium. Since the integrals in (j4.8|) have to be taken over all space 
in that case, the integration variables r" and r'" can coincide with Va, so that the 
validity of the localized form (j4.14|) of the effective susceptibility is not guaranteed. 
This comes as no surprise, as local-field effects have to be taken into account in such 
a situation. We have to conclude that (|4.15|) cannot be used as such to determine the 
average decay rate of an excited atom in an infinite medium of absorptive spheres. 
However, one is often interested in a scattering medium of finite extent, in which 
the centers of the spheres are confined to a volume V (with V/a^ ^ 1), while the 
excited atom is situated outside V. This configuration will be considered in the next 
section. The localized effective susceptibility is a useful concept in that case and 
the expression (|4.15p for the average Green function can be employed, albeit after a 
suitable modification. 



5. Finite media and surface effects 



For a finite medium the expression for the average Green function has to be modified 
so as to include finite- volume effects. We start again from (14. 2p and take the average 
over the positions of the centers of the spheres, which must be inside V. In lowest 
order of the susceptibility one encounters the average (/(r)), which for a uniformly 
distributed set of spheres with centers in V is given by 



(/(r)> 



dr' 9{a — |r — r'|) — nuo^y(r) + n 



dr' 6{a - 



'|)-«o^y(r) 



(5.1) 

instead of (|4.3|) . The step function Ov{v) equals 1 for r inside V and vanishes elsewhere. 
The expression between square brackets differs from only for positions r that are 
close to the surface of V , at a distance less than a. Assuming the surface to be 
approximately fiat on that scale, one may write r as -I- /i n, with a position vector 
at the surface and n a unit vector normal to the surface at and pointing outwards. 
In that notation one finds 

(/(r)) ^nvQ 6ly (r) +n[\vQ e{h) - tt h + \ ti h^] 6{a~\h\) (5.2) 

with e{x) = 0{x) — 9{—x). As a consequence the contribution of (G) that is linear in 
the susceptibility x{z) gets the form 

-^nv^xiz) / dr"Go(r,r",z).Go(r",r',z) 



■nx{z) I dS" I dh" [^voe{h")-Tra^h" + \Trh 



//31 



xGoir,r': + h" n" , z) ■ Go{r': 



,z). 



(5.3) 



The first term is the bulk contribution. It has the same form as the term linearly 
dependent on xiz) in (I4.15p . with the integration extended over V only. The second 
term is the surface contribution. Here dS" is a surface element at r", with a local 
normal unit vector n". 
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The surface contribution in (j5.3p may be evaluated as follows. Let us assume 
that both r and r' are far outside the volume, so that both |r — r"| and |r' — r"| are 
much larger than the wavelength (which itself is much larger than the radius of the 
spheres) . In that case the dependence of the Green functions Gq on h" is determined 
by a phase factor, as follows from (|4.11[) . As a consequence one may write 

Go(r, rl' + h" n", z) ■ Go(r'; + n", r', z) = 

= Go(r, r;,', z) ■ Go(r;,', r', z) e"'- (5.4) 

with Ss and e'^ being unit vectors in the direction r — r" and r' — r", respectively. 
For z — w + iO and q = uja/c ^ 1, as before, the exponential can be expanded. 
Subsequently, upon evaluating the integral over h" in (|5.3[) we arrive at a surface 
contribution of the form 

,2 rS 

■nvox{z)aq d^" n" • (e, + el.) Go(r, r^', z) • Go(rl', r', z) . (5.5) 

With the use of Gauss's theorem, the surface integral can be transformed to 
a volume integral. Since r and r' are both far from the surface, the ensuing 
differentiation operator can be taken to act on the phase factors in the Green functions 
only. Carrying out these differentiations we finally arrive at the following expression 
for the surface contribution in (G) that is linear in the susceptibility: 

-^nvoxiz)q^ J dr"(l + e.e')Go(r,r",z).Go(r",r',z) (5.6) 

with e being a unit vector in the direction of r — r" and an analogous unit vector 
e'. The integrand vanishes for all points r" that lie on the line connecting r and r', 
i.e. for forward scattering at the spheres, since e — — e' in that case. For backward 
scattering the correction does not vanish. 

In second order of the susceptibility one needs an expression for (/(r) /(r')) which 
takes account of finite-volume effects. Analogously to ()5.1|) we write 



z 

10c 



(/(r) /(r')) - (/(r)) (/(r')) = n j'^ dr" 0{a - |r - r"\)0{a - |r' - r"|) = 
= n0v{r) J dr"e{a-\r-r"\)6{a~\r' ~r"\) 

J dr"e{a-\r^r"\)e{a-\r'-r"\) 

-9v{r) J dr"9{a-\r-r"\)e{a-\r' ~r"\) (5.7) 

where correlation effects have been omitted, as before. The first term at the right-hand 
side leads to a volume contribution. After proper localization one finds an expression 
of the same form as the term of order [x(z)]^ in ()4.15|) with (|4.14|) . with the integration 
extended over V. 

The remaining terms at the right-hand side of (|5.7p vanish when r and/or r' are 
far from the surface. In fact, one may rewrite 6'y(r") — Ov{r) as [1 — 0y(r)]6'y (r") — 
^V'(r) [1 — 0y(r")], so that r and r" must be on different sides of the surface. Since the 
0- functions in (|5.7|) imply that these positions can at most be a distance a apart, they 
are within a distance a from the surface. As a consequence, the contribution of the 
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second term in (j5.7l) to the average Green function is a surface term. After a suitable 
change of variables it can be written as 

^ n [x{z)]^ dS" j dh" dS'" J dh'" Go(r, r^' + h" n", z) ■ 

•Go(r;' + h" n", rl" + h'" n'", z) ■ Go(r';' + h'" n'" , r', z) F(r';, v'^', h" , h'") . (5.8) 
The function F is defined as 

F{r,,r'^,h,h') ^ J ds" j dh" [~e{-'h)e{h" + \{h + h')) + e{h)e{-h" -\{h + h'))\ 

X 6* (a - |r, + i (/i - /i') n - r'; - /i" n|) 61 (a - Ir^, - i (/i - /i') n - r^' - /i" n|) . (5.9) 

The normal unit vectors at r^, r'^ and r" can be taken identical since these positions 
are at most a distance 2a apart. For the same reason the second Green function in 
(|5.8p can be replaced by its short-range approximation: 

Go(r, 0, z) ^ ^ V^, (l - 3^) + ^ ^(r) I (5.10) 

as follows from ()4.1ip . After substitution of this expression and of ()5.9p in (|5.8p the 
contribution of the delta function in (|5.10p can be evaluated along the same lines as 
before. One finds the following on a par with (|5.5|) : 



- i ^ « ^'o [x{z)? a q I' dS" n" • (e, + e'J Go(r, r^', z) ■ Go(r';, r', z) . (5.11) 

The evaluation of the contribution from the dyadic part of ()5.10|) is more 
complicated. Some of the details are given in appendix A. The result is 

- i ^ « «o [X{z)f a q I' dS" Go(r, r^', z) • (-§ I n" • e', + e', n") • Go(r;', r', z) . 

(5.12) 

A further term of second order in x(^) arises from the uncorrelated part 
(/(r)) (/(r')) of (/(r)/(r')), as given in (|5.7p . Since it is proportional to it is 
negligible for a dilute set of scatterers. 

The complete set of terms that result from surface effects in second order of the 
susceptibility x(z) is found by adding (|5.11|) and (|5.12|) : 

-i^^fo [X{z)?aq I' dS"Go{r,r':,z)-[l\n" ■{5e,+e',)+e',n"] ■ 

•Go(r';,r',z). (5.13) 
As before we may use Gauss's theorem to write this expression as a volume integral: 

,2 pV 



z 



25^2 " "0 1^ J dv" Go(r, r", z) • [I (1 + e • e') + e' e] • Go(r", r', z) . (5.14) 

As in (j5.6p the integrand vanishes for forward scattering, since e = — e' in that case 
and Go(r, r", z) • e = for large |r — r"|. 

In conclusion, we have found an expression for the average Green function 
(G(r,r',z)) of a dilute set of spherical scatterers inside a volume V. The expression 
is valid up to first order in the density and second order in the susceptibility and 
for positions r and r' far outside V. Its bulk part is given by (|4.15p (with volume 
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integrations extended over V), while the contribution from the surface is the sum of 
((^ (or dnSl)) and (fUJ^ (or ^A^). The complete resufr is 

(G(r,r',z)) = Go(r,r',z) 



■nvoxiz. 



J dv" [1 - i q2 (1 + e • e')] Go(r, r", z) • Go(r", r', z) 



+ ^nvo[x{z)?l dr"Go(r,r",z) 
•Go(r",r',z). 



I ( 2g.g, 

3 3 25 



2i 



25 



e'e 



(5.15) 

To check this expression we use it to derive the electric field generated by a 
source far away from the volume V. The Fourier component E(r,a;) of the electric 
field follows from the current density component J(r, w) of the source as 



E(r, w) = -i ^0 w J dr'G(r, r', w + iO) • J(r', oj) . 



(5.16) 



If the source is such that in the absence of the medium the electric field is a plane 
wave with Fourier component Ei(r, w) = Eq e^- exp(ik-r), with being a polarization 
vector, the average of the full electric-field component E{r,uj), including the response 
of the medium, is obtained from (|5.15l) as 

(E(r,w)) =E,(r,cj) 

, ,2 nV 



■nwox(w + iO) 



dr" 



l-^q^l-e-k) Go(r, r", u; + iO) ■ E,(r", uj) 



^2 fV 

+ —nvo[xi^ + m^ J dr"Go(r,r",c 



, I 1 1 2 1 2 

g H q I 

'3 3 ^ 25 ^ 



+ iO) • 
f 2i 3 



25^ 



E.(r", 



(5.17) 



As before, e is a unit vector in the direction r — r". The unit vector e' in (|5.15p 
could be replaced by minus the unit vector k in the direction of the wave vector of 
the incoming wave. The expression found here is consistent with that obtained for 
the average scattered field from a collection of dielectric spheres in Mie theory, as is 
shown in appendix B. 

The above derivation of the average Green function in the presence of a finite 
volume filled with absorbing scatterers clearly shows how in general both bulk and 
surface effects contribute in producing the complete result. A naive treatment in 
which the surface effects are neglected does not yield the correct answer, when the 
spheres are of a finite extent. The surface contributions account for the coarseness of 
the surface, which arises from the fact that some of the spheres may protrude. These 
protrusions give rise to specific terms in the average Green function that do not occur 
for an infinite medium. 



6. Atomic decay near a half-space of absorptive scatterers 

We consider an excited atom in the presence of a medium of absorptive scatterers that 
fills the complete half-space z < 0. The atomic position is (0, 0, Za), with Za > 0. We 
assume ZaOJa/c ^ 1, so that the results of the previous sections can be applied. In 
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the electric-dipole approximation the average decay rate fohows from l|3.10p by taking 
the average over the position of the scatterers: 



r = - 



• Im(G(r + iO)) • (ff|M|e) . 



(6.1) 



At the right-hand side we substitute expression (j5.15p for the average Green function. 
The leading term yields the standard vacuum decay rate 

, ,3 



To 



\{e\t^\9W 



(6.2) 



The next term in (|5.15l) leads to a first correction in {G{ra,ra,uja + iO)) of the 



for 



■nvox 



(cja + iO)(l- fg^) / drGo(ra,r,cja + iO) •Go(r,r„,a;a + iO) 



(6.3) 



since e = e' in the present case. The integral is a diagonal tensor, with equal xx- and 
yy-components, and a zz-component that is different. For the xx- and yy-components 
we find the following upon substituting the long-range form of the vacuum Green 
function (I4.1ip and using cylinder coordinates: 



dz 



dpp 



2z^ + 

167r(z2 +p2)2 



exp 



2i 



2N1/2 



(6.4) 



Introducing the new variable t = [z"^ + P^Y^"^ / Za instead of p and carrying out the 
integrals, we get 

^ [^^E,{u)~E,{u)-\E,{u)] (6.5) 
with w = — 2iza {iJa + iO)/c and with the functions 

/oo -xt 
dt— (6.6) 

for a non-negative integer n and for x in the right half of the complex plane. Since for 
large these functions are given by their asymptotic expansions 

n n(n + 1) 

-+ 2 +■ 

X X^ 



Enix) 



1 



(6.7) 



the XX- and yy-components of the second term in (|5.15p read 



1 / 2 

— n vq xi^^a + iO) 1 - - 
327r V 5 



(6.8) 



A similar calculation leads to the conclusion that the corresponding zz-component is 
inversely proportional to z^ so that it is small for large values of z^ uja/c. 

The last term of (|5.15|) can be evaluated in an analogous way. Upon using the 
identity Go(r, r", oja + iO) ■ e = for large |r — r"| uJa/c, the xx- and yy-components 
are found as 



1 r / •nM2 28 2 2i 



(6.9) 



while the zz-component is inversely proportional to z^, as before. 
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Collecting all results, we have found the following expression for the average decay 
rate of an excited atom in the presence of a semi-infinite medium of absorbing spherical 
scatterers: 

3 



(r)=ro 



xlm 



1 - 



16 
2 



nvo ro,_ 

1 
3 



28 



2i 



75 
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Xi^a + iO) 



(6.10) 



for large values of the dimensionless distance — ZaCOa/c between the atom and the 
medium. Here Tq_± is given by ()6.2|) . with fj, replaced by the projection fij^ of fi on 
the plane parallel to the interface of the medium. 

In deriving the above result we have taken due account of the surface effects. 
The terms of order would have been determined incorrectly, if the contributions of 
section [5] had been missed. As we have seen, the latter contributions correct for the 
fact that protruding spheres cause surface coarseness. 

Expression ()6.10p is the main result of this paper. It shows the interplay of 
absorption and scattering effects in the modification of the average decay rate. For 
spheres that do not absorb at the atomic frequency, so that xi'-^a) is real, the average 
decay rate can be rewritten as 



(r) = To - —nvoTo_± 
Id 



28 



Xi^a) 



sin(2Ca) 

Ca 



COs(2Ca) 
Ca 



.11) 



On the other hand, if absorption plays a role, whereas scattering effects can be 
neglected (as is the case for spheres with g w 0), one has 



(r)=ro 



16 



nvo To,. 



+ 1 



1. 



Xr(Wa) - 7;[Xri^^a)] 



sin(2Ca) 



COs(2Ca) 
Ca 



3.12) 



with Xr(w) and Xi{'^) being the real and imaginary part of x(w-l-iO), respectively. The 
decay rate for an excited atom in front of an absorbing dielectric half-space without 
scattering has been determined before, as we noted in section From the results in 
[171 [T^ one derives (after a few minor amendments) a decay rate for large Ca that 
coincides with (j6.12p . when the effective susceptibility (|4.14p for q — >■ is introduced. 

In figure [T] the decay rate correction function, defined as /(Ca) — ~16((r) — 
ro)/(3n?;oro.i^), is given as a function of Ca for two specific choices of the parameters 
q and xi^a), corresponding to a purely scattering case, with a decay rate given by 
(|6.1ip . and to a purely absorbing case, with decay rate (|6.12p . Both curves show a 
characteristic interference pattern. Depending on the precise location of the atom, 
the decay rate is either enhanced or reduced with respect to the vacuum decay rate. 
The effect is larger for absorbing spheres than for scattering spheres. Moreover, the 
phase of the damped oscillations of the two curves is different. For absorbing spheres 
the positions of the extrema are somewhat nearer to the interface than for scattering 
spheres. 

Interference fringes in emission processes due to refiection at an ideal mirror have 
been observed experimentally j29j . The modification of radiative properties near a 
dielectric medium has been determined experimentally a few years ago as well |22) . 
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Figure 1. Decay rate correction function /(Ca) = — 16((r) — ro)/(3nt;o Tq^x) 

for a medium with scattering spheres (with q = 0.5, x{<^a) = 0.5, ) and for a 

medium with absorbing spheres (with q = 0, xi'^a) = 0.5 + iO.5, — • — ). 



It would be interesting to mcasiire the influence of absorption and scattering in the 
medium on the emission processes and to compare the results with those found here. 

7. Conclusion and outlook 

In this paper we have shown how both scattering and absorption effects can play a role 
in the decay of an excited atom in the vicinity of a dielectric medium. Since atomic 
decay is essentially a non-classical phenomenon, a consistent treatment requires the 
use of a quantum-mechanical description for the scattering and absorbing dielectric 
medium, for the atom and for the electromagnetic fields through which they interact. 
As we have seen, a convenient model that suits these requirements is furnished by 
an inhomogeneous damped-polariton model for a set of absorptive dielectric spheres 
that scatter incoming light. In contrast to what one might expect, the average 
bulk properties of such a granular medium are in general not sufficient to account 
for all effects of scattering on the atomic decay, at least for spheres with a finite 
size as compared to the atomic wavelength. In fact, under these circumstances the 
effective dielectric constant for the bulk does not yield all information on the scattering 
processes. A subtle surface effect in the scattering contributes to the change of the 
decay as well. Only when this surface contribution is taken into account does one 
obtain the complete expression for the modified decay rate. 

In our treatment we have confined ourselves to a description of a dilute medium 
in which multiple-scattering effects are negligible. Moreover, we have considered only 
the first few terms in an expansion of the scattering amplitudes with respect to the 
ratio of the spherical diameter and the atomic wavelength. It would be interesting to 
see whether the above findings about the importance of surface effects hold as well 
when the medium gets denser or the spheres bigger. As a further simplification of 
our discussion we have assumed that the distance between the atom and the granular 
medium is large as compared to the wavelength, so that only the leading term in a 
long-range expansion of the decay rate had to be retained. For smaller distances the 
analysis of the surface effects gets more complicated. 
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Appendix A. Surface effects in second order of the susceptibility 

The surface contribution to the average Green function in second order of the 
susceptibility has been written in ()5.8|) . It contains a product of three vacuum Green 
functions, the second of which is given by its short-range approximation (jS.lOl) . The 
contribution from the delta function in this Green function has been determined in 
the main text. In this appendix we shall show how the contribution from the dyadic 
part in (|5.10p can be evaluated. 

Substitution of the dyadic part of (|5.10p in (|5.8|) with (|5.9I) yields the following 
expression: 

nixizW'P j dh" j\s"' j dh"' Go{r, r'^ + h" n",z)- 



47rc2 



. / 1 - \(r" - r'") (r" - r'") + (h" - h'")^ n" n' 

+ {h" h'") {r': - r^) n" + (h" h'") n" (r',' - r")] 



.Go(r';' + h'" n", r', z) F(r';, r^, h" , h'") . (A.l) 

The surface elements dS" and dS'" are located at r" and r"', respectively. The 
normal unit vectors at these two positions are almost equal and have been denoted 
by n". The principal value sign indicates the exclusion of a small sphere around 
r" + h" n" in the integrations over r"' and h'" . Furthermore, F{ys, r^, /i, h') stands for 
F(rs, r^, /i, h')/[\vs — r'^.^ + {h — h'Y]^/"^, while the variable z equals w + iO, as before. 

As in (|5.4I) , the two Green functions in (jA.l[) can be expanded (in their second or 
first argument, respectively) around their values at r", if both |r — r"| and |r' — r"| 
are large compared to the wavelength, as we have assumed before. Up to first order 
in wa/c, the ensuing phase factor can be expanded as 

1 - i - /i" n" • e, - i - h'" n" • e' + i - (r" - r'") • e' (A.2) 

c c c 

since h" , h'" and |r" — r"'| are all of order a at most. As before, e^ and e'^ are unit 
vectors in the direction of r — r" and r' — r", respectively. The product of (|A.2[) and the 
expression between curly brackets in (|A.ip contains all information on the dependence 
of the integrand on r" — r"'. This product may be replaced by the effectively equivalent 
form 

; |r'; - r'/'p + {h" - h'"f] (l - i ^ h" n" ■ e. - i ^ h'" n" • e',) (1-3 n"n") 
_| i f (/," _ h'") rr P (e; n" + n" e^ - 2 n"n" n" • e^)} 

x[\r':^r':'\' + {h"-h"')T' (A.3) 

up to first order in wa/c, when use is made of the rotation symmetry of the integration 
over r'" in the planar surface through r" and orthogonal to n". This symmetry implies 
that the product (r'/ - r^,") (r'/ - r'/') may be replaced by i |r;,' - r;,"^ (I - n"n"). 

To simpliiy the multiple integral in (jA.l[) we introduce, instead of h" and h'" , 
their sum and difference as new integration variables. Subsequently, we carry out 
the integral over h" + h"', at fixed h" — h'" . Inspecting the result, one finds that in 
leading order of uja/c the expression (jA.l|) vanishes owing to the odd parity of the 
integrand in its variable h" — h'" . The contribution from the next order in uja/c does 
not vanish. It can be determined by first considering the integrations that are hidden 
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in the definition of the function F, as given in (|5.9p . These lead to the following 
integrals: 

Iiirs,h) = f dS' j dh'e[a^ -\\v,-r'f -{\h-h'f] 

X0[a2 _ 1 1 + rlf - (i /i + h' f] (A.4) 

J2(r„ h) = dS' j dh' h'^ 9[a^ - |i - r;|2 - (i /i - h'f] 

x^?[a2_|lr,+r^|2-(l/^ + /^')2] (A.5) 

with ft, > 0. In writing these integrals we have chosen the origin of the coordinate 
system to be situated at the surface. Both of the integrals vanish for r = [jr^p + 
/i2]i/2 > 2 a. 

In terms of the above integrals the contribution (jA.l[) to the average Green 
function becomes in leading order of uja/c: 

- i^ n [x{z)r V f dS" J' dS'" J dh'" [|r'; - r^'f + {h" ft'")^]-^/^ Go(r, r^', z) ■ 

• [[-i |r;' - r'^f + [h" h"'f] {[\h - \{h" - h"'f h] n" ■ (e, + e^) 
+ 1 _ „" . (e^ _ el)} (I - 3 n"n") 

+ 1 {h" - h"r W: ~ r';f h (e; n" + n" e', - 2 n" n" n" • e'J] • Go(r';, r', z) (A.6) 

with Ii and I2 depending on r" — r"' and h" — h'" . 

In order to proceed we have to determine explicit expressions for the integrals 
Ji and 12- The first is the standard overlap integral, which has been encountered in 
(|4.4p . It equals vq c(r), with c(r) given in (|4.7p . The integral I2 is the second moment 
of the overlap integral. Choosing cartesian coordinates in such a way that the normal 
to the surface at the origin points in the direction of the positive z-axis, and that 
equals (a;, 0, 0) with a; > 0, we may write (|A.5[) as 

dx' dy' dh'h'H[a'^~-{^x-x'f-y''^-{^h-h'f] 

-00 J —00 J —00 

•K6[a^- {\x + x'f - y'2 - (ift + h'f] . (A.7) 
The integral over y' is trivial, with the result 

dh' h'^ / dx' [a^ - (ix + x'f - (i/i + h'f] 

-00 J —hh' I X 

xe[a^- {\x + x'f - {\h + h'f] . (A.8) 

The 0- function constrains the integrations over x' and h' . This constraint 
interferes with the bounds on the integration written explicitly in (jA.Sp . By a 
geometrical analysis one finds that the interference depends on the sign of the 
combination x^ + h^ — 2ah. In fact, for 2 a ft, < x^ + h^ < 4 the upper limit of the x'- 
integration is effectively finite, while the lower limit is left unchanged. Furthermore, 
the domain of the ft'-integral is found to be constrained to values |ft'| < M, with 
M = X [(a^ — jx"^ — J h^)/{x^ + ft2)]^/2. Hence, after a shift of the x'-variable one has 
the following for 2 a ft < a;^ + ft^ < 4 a^: 

/2(r„ft)=4 / dft'ft'2 / dx' [a^ -x'^ ~{\h + h'f]' . (A.9) 

J-M J -hh' /x+x/2 
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a - T a; 



1 ™2 1 7 2 



1/2 



-i [a^ -{\h + h'f] arcsin 



hh' 



1/2 



(A.IO) 



After a partial integration in order to get rid of the arcsine function, the integral over 
h' can be carried out as well. The final result is 

J.(r., h) = K (-1 - 5 Z^') + - ^ - 1^ /^^) 

-3M^'-2io^''^'+9M^''^'+4io'^']+l^'^«' (4a2 + 5/.2) (A.ll) 
for 2a/i < x2 + < Aa^. 

In the other case x^ + h? < 2 ah the double integral is the sum of two contributions 
with different bounds: 



dh' h'^ / dx' [a^ - x'^ - {\h + h'f] 

-M J -hh' /x+x/2 



-h/2+a 



dh' h'- 



M 



/x+x/2 
[a^-(h/2+h'f]^''' 



dx' [a^ -x''^ - i^h + h') 



21 1/2 



(A.12) 



-[a-^-(h/2+h'Y\^>^ 

Upon evaluating the x'- and the ft'-integral we arrive at a result that is found to be 
identical to that given in (jA.llI) . 

Employing spherical coordinates in (|A.lip . with h ~ r cos0 and x = r sin0, we 
may write the second moment of the overlap integral for both cases as 



l2{rs,h) = voa"^ 



3r 
16^ 



(1 + cos^ 



COS 



1280a5 



(3-5 cos^ 



4a2 
e{2a-r) . 



(1 — 3 cos^ ( 



(A.13) 



Having obtained explicit expressions for Ii and /2, we return to (|A.6|) . 
The integrals over r"' (with surface element dS'") and h'" can be calculated 
straightforwardly upon introducing spherical coordinates and performing the angular 
integration first. In this way we arrive at the result 



25c' 



■nvo [x{z)?aq j' dS" Go(r,r;', z) .(-fin", e', + e^, n" + n" e'J 



•Go(r';,r',z) 



(A.14) 



This is equivalent to (|5.12p . since one may use the identity e^, • G(r",r',z) = for 
distances |r' — r"| that are large compared to the wavelength. 



Appendix B. Scattering fields in Mie theory 

Electromagnetic scattering from a dielectric sphere has first been treated by Mie [30] 
and reviewed subsequently by several authors |31l I32j . If a linearly polarized incoming 
plane wave, with the wave vector k in the direction of the positive z-axis, with the 
polarization vector e„ along the x-axis and with the amplitude i?o, impinges on a 
dielectric sphere with a radius a, with centre at the origin and with a dielectric constant 
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e = 1 + Xi the components of the scattered electric field in the far-field region have 
the form 



El;>{r,e,^) ^Eo— cos(^ V(-i)'^ [BlTe{cos9) + BY'ne{cosl 

rZT 



[r, if) = ~Eo — sin ip ^^(-i)' [5f ^^cos 9) + n{cos 0)] 



(B.l) 



with spherical coordinates r, 9 and (p. The angular functions are defined in terms of 
associated Legendre polynomials as 



TTg (cos 9) 



1 



smf 



(cos 6*) , (cos 6*) = — P/ (cos t 
dO 



The electric and magnetic multipole amplitudes read 
„p ^ .,+1 2£+l Nl 

with p — e, m. The numerators and denominators are given as 

N! = e [{l+l)j,(q) - qn+i{q)] Mq') - [{i+l)ji{q') - q' Je+l{q')] Mq) 

NT = [(^ + ^)3t{q) - q]i+i{q)] Mq') - [(^ + i)Mq') - q' je+i{q')] Mq) 



(B.2) 
(B.3) 



{£ + !) hf\q) -qh^llM] Mq') - [{t + i)Mq') - q' u+M)] hf\q) 

i + \)hf \q)-qhfl{q)\ Mq') - [{^ + ^) Mq') - q' u+M)] hf \q) (b.4) 



with spherical Bessel and Hankel functions depending on q = ka and q' = -Je q. 
For small values of q, the first few multipole amplitudes get the form 



3 + X 



5 3+x 3 3+x 



1 5 X 

■q 



18" 5 + 2x 



1 30*^ ^ 



(B.5) 



up to the order q^. When x is small as well, the first two of these can be written as 



Bl^-q^X 



1 1 1 2 

1 — X — q 
3^ 5^ 



2i 



q^x 



90 



9'^ X 1 



(B.6) 



up to the order x^. Upon substitution in (jB.ip the far fields are found as 



E^"^ {r,9,ip) ^ Eo^q-^XCOSip [l-^q^] cos + -f- q^ cos 



+X 



kr 

* ( 1 - 2i 



15 



cos 61 - — 9^ (2 cos^ 61-1) 
75 



e f 1 / 1 \ 1 

4^)(r,0,^) = -i;o — g^^^sin,^ |- h +_^2 ^Qg^ 



1 



2i 



-n ( 1 -9 - -9 ) - cos 



1 



(B.7) 



Atomic decay 



20 



up to the order and x^- In vectorial notation this expression may be rewritten by 
introducing the long-range form of the vacuum Green function (14.11^ : 

,2 



E 



+X 



c2 
1 1 



2i 
9"' 



iO) 

1 



1 9 



I f • k + kf 



i,^f.k 



E,(0) 



(B., 



with the spherical volume vq = 4 tt a'^/S and with the unit vectors r = r/r and k = k/fc. 
This expression for the scattered electric field is consistent with that found in (|5.17p 
for the average field due to scattering from a set of Mie spheres. 
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